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Abstract. - Phase transitions induced by varying the strength of disorder in the large-g 
state Potts model in 3d are studied by analytical and numerical methods. By switching on the 
disorder the transition stays of first order, but different thermodynamical quantities display 
essential singularities. Only for strong enough disorder the transition will be soften into a 
second-order one, in which case the ordered phase becomes non-homogeneous at large scales, 
while the non-correlated sites percolate the sample. In the critical regime the critical exponents 
are found universal: fj/i> = 0.60(2) and v = 0.73(1). 



Real materials are generally three-dimensional and contain some sort of quenched disorder, 
whose effect on the properties of the phase transition is an intensively studied question. If 
the phase transition in the non-random system is of second order, relevance or irrelevance of 
disorder can be decided due to the Harris criterion [1] and also perturbative field-theoretical 
calculations are possible [2]. For first-order transitions in the pure system, however, we have 
only a more limited knowledge. In 2d any amount of continuous disorder is sufficient to soften 
the transition into a second-order one [3,4], which is studied in detail numerically [5-7]. In 
3d, however, there is no such a rigorous criterium. Here experimental work [8] and numerical 
simulations [9-11] show that for strong enough disorder the phase transition becomes second 
order. In the numerical work mostly the q = 3— and 4-state Potts models were studied by 
Monte-Carlo simulations, both by site [9,10] and bond dilutions [11]. 

In the present work we revisit the problem of the effect of disorder at first-order phase 
transitions in 3d and address mainly such questions which have not been treated previously. 
We consider the disorder as realized by random ferromagnetic bonds, which contains as a 
limiting case the bond dilution. The problem we investigate is the Potts model in the large-g 
limit, which allows us to make partially analytical calculations and in this way to shed light on 
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the physical mechanism responsible of the softening of the transition. We have also performed 
extensive numerical calculations, in which exact results on large finite samples are averaged 
over disorder. In particular, we addressed the question of universality of the second order 
transition and the location of the tricritical point. 

We start to define the Hamiltonian of the g-state Potts model [12]: 

H = Jij^i,<Tj) (1) 

(id) 

in terms of the Potts-spin variables, <7j = 0, 1, • • • , q — 1, at site i and the summation runs 
over nearest neighbors of a simple cubic lattice. The ferromagnetic couplings, Jjj > 0, can 
take two values, J\ = J(l + 5) and J 2 = J(l — 5) with equal probability [13]. The strength of 
disorder is 8 and the dilute model is recovered for S = 1. As usual for the random model we 
calculate disorder averages of the free energy and the correlation function. 

For a given realization of disorder the partition function is given in the random cluster 
representation [14]: 

z =y,i c(g) n [e^'-i] ( 2 ) 

G ijeG 

where the sum runs over all subset of bonds, G, c(G) stands for the number of connected 
components of G and we use the reduced temperature [7], T — > Ting = 0(1) and its inverse 
as (3 — > p/lnq. In the large-*/ limit, where q@ Jij ^> 1, the partition function can be written as 

z=y j ^ (g \ m = c(G)+pj2 J v ( 3 ) 

GCE ijeG 

which is dominated by the largest term, <p* — maxg <fr(G) , and the degeneracy of the optimal 
set G* is likely to be one. Thus the study of the large-g state Potts model is reduced to 
an optimization problem, whose solution (for a given realization of disorder) depends on the 
temperature. The free-energy per site is proportional to 4>* and given by — (3f = <fi* /N, where 
N stands for the number of sites of the lattice. From the point of view of magnetization and 
the correlation function the largest connected cluster of the optimal set, C, plays a crucial 
role. The magnetization, m, is given by the fraction of sites in C, and its disorder average is 
m > in the ordered phase, T < T c . In the disordered phase, T > T c , C has a finite linear 
extent, £ < co, which is proportional to the average correlation length in the system. At the 
phase-transition point, T = T c , if the transition is of second order, C is a fractal and its fractal 
dimension, df, is related to the magnetization scaling dimension, x m , as d = df + x m . Here 
be expressed by the standard critical exponents as x m — [3 /v. On the contrary, if the 
transition is of first order there is phase coexistence at T — T c and the correlation length is 
finite. This type of formalism has been used by us previously in 2d [7], and we refer to this 
work for a detailed presentation of the method. 

The solution of the problem is the simplest for the non-random model, 5 = 0, when there 
are only two homogeneous optimal sets. For T < T c (0) it is the fully connected diagram with 
a free-energy: —(3Nf = 1 + N/3Jz (z = d = 3 is the connectivity of the lattice) and for 
T > T c (0) it is the empty diagram with —0Nf = N. Consequently the transition point is 
located at: T c (0) = Jz/{1 - l/N) and the latent heat is Ae/T c (0) = 1. 

Introducing disorder into the system, 5 > 0, new type of non-homogeneous optimal di- 
agrams will appear close to T c (0). Let us consider first T > T c (0), where the new optimal 
set contains connected clusters of linear size, I, embedded into the empty diagram. Here we 
refer to the treatment of this problem in 2d [7] when close to T — T c (0) the excess quantities 
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over the homogeneous optimal sets are i) a bulk term due to disorder fluctuations, A/;,, and 
ii) an interface term due to missing bonds, A/ s ~ f3JS, where the size of the interface is 
S ~ l d ~ x . For Gaussian fluctuations A/^ 5 ~ 8j3 JV 1 / 2 , where V ~ l d is the volume of the 
cluster. Consequently for d < 2 normal fluctuations of disorder will lead to the creation of 
clusters of arbitrary size and for any weak disorder there is no phase coexistence and the 
transition is of second order [3] . 

On the contrary for d > 2, clusters are created only due to extreme fluctuations, when 
say the cluster consists of only strong bonds, thus A/^~ ~ (3J5(V — cS), where c = 0(1) is 
a geometrical factor. The existence of such large cluster is exponentially rare, its density is 
given by — In p + ~ V — cS. For a small 8 only large clusters can be formed and we obtain 
the possible size from the condition, A/^~ > A/ s , as I > l+(8) ~ 1/8 [16]. Finally, the free 
energy of the non-homogeneous optimal set relative to the empty graph is given by the sum 
of contribution of the possible clusters: 

-PNf+~N + N J2 ( A /6 + (0 - A/ s (/))p+(0 , (4) 

l>l + (S) 

and it is dominated by the contribution with I = l+(8). 

A similar analysis can be performed for T < T c (0), when the new inhomogcncous optimal 
set is obtained from the fully connected graph by creating clusters of isolated sites. In this case 
a cluster of size, I, has V + cS weak bonds, a bulk gain of A/ b ~ ~ f3J8(V + cS) and appears 
with a density, — In p_ ~ V + cS. Note that there is a difference in the sign of the surface 
term comparing with the high-temperature case. The limiting size of clusters, I— (8) can be 
estimated as above and the free energy of the optimal set can be written in an analogous form 
with eq. (QJ: 

-f3Nf_~l + N(3Jz + N (A/ fe -(0-A/ s (0)p-(0, (5) 

l>l-(5) 

and it is dominated by the contribution with I = l-{8). The asymmetry in the densities in the 
two phases, p+ and />_, leads to a shift in the critical temperature, ln(T c (<5) — T c (0)) ~ — <5 -3 , 
and similarly to a reduction of the latent heat and the jump of the magnetization at the 
transition point: 

ln(l- Ae/T c ) ~ ln(l - Am) ~ ~ . (6) 

Thus the phase transition stays first order and as the disorder is switched on there is an 
essential singularity in the thermodynamical quantities as a function of 8~ 3 . 

For finite values of disorder the optimal set (and thus the free energy and the magneti- 
zation) is exactly calculated by a combinatorial optimization algorithm [15], which works in 
strongly polynomial time. In the numerical calculations we studied periodic systems of cubic 
shape with a linear size L = 16,24,32 and in some cases L — 40. For the largest size the 
computation of the optimal set for a single sample typically took from 5 to 20 hour CPU time 
in a 2.8GHz processor. (Since the breaking-up length [16] at 8 = 0.40 is Lt, r « 40 we restricted 
ourselves to 0.40 < 8 < 1.) The number of realizations were several thousands, for the largest 
size several hundreds. 

We start to show in fig. ^ the numerically calculated phase diagram as a function of the 
strength of disorder. Here the phase boundary between the ordered and disordered phases 
is obtained from the condition that the largest cluster in the optimal set starts to percolate 
and the finite-size results are then extrapolated, see the inset of fig. Q The transition points 
obtained by this method agree, within the error of the calculation, with those calculated from 
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Fig. 1 - Phase diagram of the random-bond Potts model for bimodal disorder. The calculated 
transition points are denoted by crosses the size of which is larger than the error of the calculation. 
In the ordered phase in the optimal set the non-connected sites are percolating for S > S pr and 
T pr (8) < T < T C (S), see eq. J7J. Numerical results in fig. ^indicate, that the first- and second-order 
transition regimes are separated by a critical disorder, S c , which corresponds to the border of the 
percolation regime, S c = S pr . Inset: Distribution of the finite-size transition points identified with the 
temperature at which the connected cluster starts to percolate. 



the position of the maxima of the specific heat. Before analyzing the properties of the phase 
transition we mention that the structure of the optimal set in the ordered phase can be of two 
different types. For weak disorder, 8 < 8 pr , or for low temperature, T < T pr , in the optimal 
set the isolated sites form finite clusters, otherwise the isolated sites are percolating. In the 
dilute model with 8=1 and J\ = 2 J the percolation temperature is at T pr (l) = J\. Indeed 
the optimal set for T < J\ contains all the strong bonds, whereas from this optimal set for 
T > J\ the dangling bonds are removed. Since the dangling bonds represent a finite fraction 
of the bonds the non-occupied sites become percolating. 

For 8 < 1 the skeleton of dangling bonds of J% couplings are decorated by weak J2 cou- 
plings, and by removing one strong dangling bond one has also to remove four weak decorated 
ones, which is possible in the temperature range: 

T c > T > T pr (8) = J(5 - 38) . (7) 

The numerical results indicate that at T pr {8) in a finite fraction of samples there is a giant 
cluster of isolated sites which spans the finite cube. 

The temperature dependence of the magnetization for different values of the disorder is 
illustrated in fig. [5] in a finite system of L — 16. For 8 < 8 C in the first-order transition 
regime the magnetization has a finite jump, Am > 0. Am(L) is estimated from the relation: 
Am(L) sa 2m[T c (L)}, and the extrapolated values are shown in the inset of fig. [3for different 
values of 8. From this we estimate the position of the tricritical disorder 8 C = 0.65(2). It 
is interesting to observe that this value, within the relatively large error coincides with the 
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Fig. 2 - Temperature dependence of the magnetization for L = 16 with different strength of disorder: 
5 = 0.5, 0.6, 0.68, 0.75, 0.8, 0.875, 0.925 from left to right. Inset: Jump of the magnetization at 
the transition point as extrapolated at 8 — 0.4, 0.5 and 0.6. In the region to the right of the arrow, 
S > 0.68, the fractal dimension is obtained df < 3 (see the method in figEJ which indicates a second 
order transition, thus Am — 0, which is shown by crosses. The line drawn over the extrapolated 
points is a guide to the eyes. 



limiting value of disorder, 5 pr = 0.658(1), where the percolation of isolated sites in the ordered 
phase starts. Physically it is very intuitive that percolation of isolated sites is the prerequisite 
of a continuous transition in the system. Indeed the correlation length in the ordered phase, 
which is measured by the linear size of the largest finite (i.e. not the giant) ordered cluster 
can be divergent, if it is embedded in an infinite cluster of isolated sites. Therefore we tend 
to conjecture that 5 C = S pr . 

In the second-order transition regime we have made the most detailed calculations at 
S = 7/8, but universality is checked at other values of the disorder at S = 0.75, 0.80, 0.925 
and 1. The fractal dimension, df, of the percolating connected cluster is estimated in such a 
way that fixing a reference point we measured the average number of points, s(r, L, T), within 
a shell around the reference point of unit width and radius, r. Close to the transition point: 
t = (T — T c )/T c -C 1 it is expected to scale as: 

s(r, L, t) = L^-^ir/L, tL 1 '") . (8) 

Since the exact value of the critical temperature is not known we can not fix the second 
argument of the scaling function as s(p,r = 0). Instead we set T = T sp (L), which is the 
temperature, where the connected cluster spans the cube of the given size and we average 
over disorder. Evidently with this choice the second argument of the scaling function, t, is 
asymptotically constant and the scaling function depends only on one parameter: s — s(r/L). 
Our scaling picture is verified in fig. [3Jin which the scaling plot of the mass in the shell is shown. 
From the optimal scaling collapse [17] we obtained for the fractal dimension, df = 2.40(2). 
The correlation length critical exponent, v, is calculated from the scaling behavior of the 
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Fig. 3 - Scaling plot of the mass of a shell of the giant cluster, for each size at the average spanning 
temperature, see text. Optimal collapse is obtained with a fractal dimension, dj = 2.40. Inset: 
Scaling plot of the magnetization as a function of the distance from the critical temperature. By 
fixing p/v — d — df = 0.60 the best collapse is obtained with v — 0.73(1) 

magnetization: m(t,L) = L Xm m(tL l t v ). From an optimal scaling collapse [17] as shown in 
the inset of fig. |21we have obtained v = 0.73(1). Note that v satisfies the rigorous bound for 
disordered systems [18]: v > 2/d. Due to strong cross-over effects we could not determine 
the tricritical exponents with sufficient accuracy. However, the tricritical correlation length 
exponent, v tr , is related [6] to the exponents of the random-field Ising model as: v tr — 
v rf 7(2 — OiRF — Prf)- Using numerical estimates [19] we obtain v tr ~ 0.69. 

In Table [I] we have collected the presently known numerical results about the critical 
exponents of random q-state Potts models in 3d, including also the Ising model (q = 2) and 
percolation, which formally corresponds to q — 1. The q-dependence of the exponents, in 
particular that of (3/v is non-monotonic. 

In our model thermal fluctuations play a negligible role and the physical properties are 
solely determined by disorder effects, which is characteristic at a so called infinite disorder 
fixed point [20]. In the two-dimensional problem an isomorphism is conjectured [7] with 
the phase transition in the strongly anisotropic random transverse-field Ising model in Id. 
However, our numerical results indicate that this type of relation does not seem to work in 

Table I - Critical exponents of the random q-state Potts model in 3d comparing with those of perco- 
lation 





perc. [21] 


<7 = 2 [22] 


q = 3 [10] 


q = 4 [11] 


q — > oo 


V 

P/u 


0.877 
0.477 


0.684(5) 
0.518(2) 


0.690(5) 
0.539(2) 


0.747(10) 
0.732(20) 


0.73(1) 
0.60(2) 
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higher dimensions. 
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